Experimental test of uncertainty relations for quantum mechanics on a circle 
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We rederive uncertainty relations for the angular position and momentum of a particle on a circle by employ- 
ing the exponential of the angle instead of the angle itself, which leads to circular variance as a natural measure 
of resolution. Intelligent states minimizing the uncertainty product under the constraint of a given uncertainty 
in angle or in angular momentum turn out to be given by Mathieu wave functions. We also discuss a num- 
ber of physically feasible approximations to these optimal states. The theory is applied to the orbital angular 
momentum of a beam of photons and verified in an experiment that employs computer-controlled spatial light 
modulators both at the state preparation and analyzing stages. 
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I. INTRODUCTION 

Despite the outstanding role that angular variables play in 
physics, their proper definition in quantum mechanics is beset 
by difficulties and requires more care than perhaps might be 
expected [ Ij, |2|, |3t] . Consider, for instance, the simple example 
of a particle moving on a circle of unit radius: the problems 
essentially arise then from the periodicity, which prevents the 
existence of a well-behaved angle operator, but not of its com- 
plex exponential, we shall denote by E. 

In quantum optics, this topic is by no means purely aca- 
demic: it turns out to be crucial for a proper understanding 
of, e.g., the orbital angular momentum (OAM) of light |i4J. 
Indeed, as put forward by Allen and coworkers j^, the 
Laguerre-Gauss modes, typical of cylindrical symmetry, carry 
a well-defined OAM per photon. Since it is surprisingly sim- 
ple to generate, control, filter, and detect OAM states of light 
experimentally, researchers have begun to appreciate their 
practical potential for classical [6, 7, 8] and quantum infor- 
mation applications ||1 E llll liJ, 111 H H . 

Intimately linked to the issue of a proper angle descrip- 
tion it is the question of the associated uncertainty relations. 
Surprisingly enough, some subtle aspects of these relations 
still remains under discussion. From previous work in this 
topic 1 16, 17, 18, 19, 20,ll|2l|2l|2l it seems clear that, if 
one insists in holding to an angle operator, special care must 
be taken when using the standard variance, since this is a non- 
periodic measure of spread that makes the angular uncertainty 
depend on the 2tt window chosen. Moreover, the associated 
commutation relation depends on the value of the angle dis- 
tribution at a point, which turns it of somewhat cumbersome 
handling. 

By precise measurements on a light beam, a detailed test 
of the uncertainty principle for angle and angular momentum 
has been recently demonstrated jzsl 126 1. The idea is to pass 
the beam through an angular aperture and measure the result- 
ing angular-momentum distribution iIztIi . In the same vein, 
we have presented experimental results f2§\ that strengthen 
the evidence that E furnishes a correct description of angu- 
lar phenomena. When a sensible periodic resolution measure 
(namely, the circular variance) is employed, the associated in- 



telligent states should minimize two inequalities (one for the 
cosine and other for the sine), and both cannot be saturated si- 
multaneously. To bypass this drawback, we have looked at the 
more physically meaningful notion of constrained intelligent 
states; that is, states that minimize the uncertainty product for 
a given spread either in angle or in angular momentum. In 
fact, they prove to be Mathieu wave functions, which have 
been attracting great interest in relation with nondiffracting 
optical fields iilM in . 

In this paper we go one step beyond and present an im- 
proved experimental setup (that uses computer-controlled spa- 
tial light modulators both at the state preparation and ana- 
lyzing stages) to verify in great detail the properties of these 
constrained intelligent states. As a byproduct, we also bring 
out that E can be associated with a feasible transformation 
(a fork-like hologram) that shifts the values of the angular 
momentum. Our formulation paves thus the way for a full 
quantum processing of vortex beams and provides a bridge 
between the classical theory of singular optics and the realm 
of quantum optics. 

The plan of this paper is as follows. In Sec. |ll]we provide 
a comprehensive quantum treatment of angular variables, in- 
cluding a discussion about the associated coherent and intel- 
ligent states, as well as various suboptimal states. A feasible 
optical realization of the system under study is provided in 
Sec.|III] with special emphasis on the detection of the angular- 
momentum spectrum. In Sec. |IV]our setup is shown and ex- 
perimental results are presented and discussed at length. Fi- 
nally, the summary of our achievements and suggestions for 
possible future upgrades are given in Sec. [V] 



II. THEORY 
A. Quantum description of rotation angles 

We consider rotations of angle (p generated by the angu- 
lar momentum along the z axis, which for the simplicity we 
shall denote henceforth as L. Classically, a point particle is 
necessarily located at a single value of the periodic coordi- 
nate (f), defined within a chosen window, e.g., [0, 2tt). The 
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corresponding quantum wave function, however, is an object 
extended around the unit circle §i and so can be directly af- 
fected by the nontrivial topology. 

If we treat as a continuous variable, the Poisson bracket 
for the angle and the angular momentum is 

{0,L}-1. (2.1) 

Direct application of the correspondence between Poisson 
brackets and commutators, suggests the commutation relation 
(in units h = 1) 

[$,L]=i. (2.2) 

One is tempted to interpret as multiplication by 0, and then 
to represent L by the differential operator 




that verifies the fundamental relation ( 12.21 ). However, the 
use of this operator may entail many pitfalls for the unwary: 
single-valuedness restricts the Hilbert space to the subspace of 
27r-periodic functions, which, among other things, rules out 
the angle coordinate as a bona fide observable Il32| [33[ [34ll . 
A possible solution, proposed by Judge and Lewis IBSll . is to 
modify the angle operator so that it corresponds to multipli- 
cation by (p plus a series of step functions that sharply change 
the angle by 2tt at appropriate points, which coincides with 
the classical Poisson bracket of L and a periodic variable. 

Many of these difficulties can be avoided by simply se- 
lecting instead angular coordinates that are both periodic 
and continuous. However, a single such quantity cannot 
uniquely specify a point on the circle because periodicity im- 
plies extrema, which excludes a one-to-one correspondence 
and hence is inc omp atible with uniqueness. Perhaps the sim- 
plest choice lf36ll37ll is to adopt two angular coordinates, such 
as, e. g., cosine and sine. In classical mechanics this is indeed 
of a good definition, while in quantum mechanics one would 
have to show that these variables, we shall denote by C and 
S to make no further assumptions about the angle itself, form 
a complete set of commuting operators. One can concisely 
condense all this information using the complex exponential 
of the angle E — C + iS, which satisfies the commutation 
relation 

[E,L]^E. (2.4) 

In mathematical terms, this defines the Lie algebra of the two- 
dimensional Euclidean group E(2) that is precisely the canon- 
ical symmetry group for the cylinder §i x M (i.e., the classical 
phase space of the system under study). 

The action of E on the angular momentum basis is 

E\m) = |m - 1) , (2.5) 

and, since the integer m runs from — oo to +oo, E isa unitary 
operator whose normalized eigenvectors are 

oo 



The intuitively expected relationship of a discrete Fourier 
pair between angle and angular momentum is an immedi- 
ate consequence of Eqs. ( 12.3b and ( 12.6b . Indeed, denoting 

^rn = {m\^i) and = {4>\^), it holds 

^ m— — oo 

(2.7) 

Jo 

There is an appealing physical interpretation beyond the 
definition of E. Whereas in the case of ( 12.2b one thinks in 
terms of complementarity between two measurable quantities, 
E primarily represents a transformation and ( 12.4b may be in- 
terpreted as the complementarity between measurement and 
transformation. On the other hand, the action of E can be cast 
also in terms of measurement, since any unitary operator may 
be generated by an appropriate Hermitian generator There 
is a twofold goal in the theory of angular momentum and its 
conjugate variable: to characterize them either as a transfor- 
mation or as a measurement. Notice that the pioneering work 
in Refs. [25] and [26] anticipated the former interpretation. 

The role of _E as a transformation is determined by the ac- 
tion on the basis states ( 12.5b . What deserves to be explained 
is the possible measurement associated with E. Although the 
vectors |0) provide an adequate description of angle, one must 
take into account that realistic measurements are always im- 
precise. In particular, the measurement of the spectrum would 
require infinite energy. In other words, the mathematical con- 
tinuum of angles will be observed always with a finite reso- 
lution. In consequence, it could be interesting to extend the 
previous formalism by including fuzzy, unsharp or noisy gen- 
eralizations of the ideal description provided hy E. To this end 
we shall use positive operator-valued measures (POVMs), that 
are a set of linear operators A(0) furnishing the correct prob- 
abilities in any measurement process through the fundamental 
postulate that llssll 

= Tr[M(0)] , (2.8) 

for any state described by the density operator g. Compati- 
bility with the properties of ordinary probability imposes the 
requirements 

A(0)>O, A(0)= At (</,), r dcp A{cf>) ^ i . (2.9) 

Jo 

In addition to these basic statistical conditions, some other 
requisites must be imposed to ensure a meaningful description 
of angle as a canonically conjugate variable with respect L. 
We adopt the same axiomatic approach developed previously 
by Leonhardt et al [39.1 for optical phase. First, we require the 
shifting property 

e"l''^ A{(j)) e-"^''^ = A(<?!> + <?!>') , (2.10) 



3 



which reflects nothing but the basic feature that an angle 
shifter is a angle-distribution shifter and imposes the follow- 
ing form for the POVM fi^l 



1 

2^ 



E 



A. 



^i{m--m') 



'\m){m'\. (2.11) 



m,m — — oo 



We must also take into account that a shift in L should not 
change the phase distribution. But a shift in L is generated by 
E since, according to (12.5b . it shifts the angular momentum 
distribution by one step. Therefore, we require as well 



(2.12) 



which, loosely speaking, is the physical translation of the fact 
that angle is complementary to angular momentum. This im- 
poses the additional constraint Am+i,m'+i — ^m.m', and this 



means that A„ 
be recast as 



In consequence, Eq. (12.11b can 



1 — — OQ 

and the conditions i2.9i are now 



Ar = A_, . 



(2.13) 



(2.14) 



Expressing E in terms of its eigenvectors, we finally arrive at 
the more general form of the POVM describing the angle vari- 
able and fulfilling the natural requirements ( 12.10b and ( 12.12b : 



A(0)= / d<j>' IC{<P')\<P + (b'){cj} + cj,' 



where 



/C(0) 



1 

2^ 



E A 

1 — — 00 



I e 



(2.15) 



(2.16) 



The convolution ( 12.151 ) shows that this POVM effectively rep- 
resents a noisy version of the usual projection 10) (f/)!, and the 
kernel /C((/)) gives the resolution provided by this POVM. 



B. Gaussian distributions on a circle 

Experience with quantum mechanics of simple systems, 
such as the free particle and harmonic oscillator, suggests that 
Gaussian states can be an important tool for a better under- 
standing of the periodic motion on a circle. Even so, it is re- 
markable that there is no clear concise definition of the Gaus- 
sian distribution on a circle and one can only find vague state- 
ments scattered through the literature. 

We do not want to enter here into a mathematical treatment, 
but rather we just try to grasp the properties that make the 
Gaussian distribution on the line to play such a key role in 
physics and that we are particularly keen on retaining when 
constructing its circular counterpart. We itemize the most rel- 
evant ones in our view: 



1 . The sum of many independent random variables tend to 
be distributed following a Gaussian distribution. 

2. All marginal and conditional densities of a Gaussian are 
again Gaussians. 

3. The Fourier transform of a Gaussian is also a Gaussian. 

4. The Gaussian distribution maximizes the Shannon en- 
tropy for a fixed value of the variance. 

The first property (subject to a few general conditions) is 
the central-limit theorem and explains the ubiquity of Gaus- 
sians in physics: the distribution of the phenomenon under 
study does not have to be Gaussian because its average will 
be. The second and third ones are responsible for the good 
properties than one assigns to Gaussian states in quantum op- 
tics. Finally, the last condition bears on the information-based 
approach to quantum theory, but strongly depends on the def- 
inition of entropy we adopt. 

In statistics there are two distributions that have been some- 
how suggested for having good properties on a circle, namely 



Pa{<P) = 



1 



■ exp[2KCOs(0 — fi)] , 



1 



OO 



fe— — oo 



exp 



1 (0-^ + 27rfc)2 



(2.17) 



R.18) 



where /„ denotes the modified Bessel function of first kind. 
The first one is known as the von Mises distribution, while the 
second is the wrapped Gaussian. By a trivial application of 
the Poisson summation formula we can express the latter as 



P(t(0) 



1 

2^' 



1^3 U - M 



1 



where 



MC\q)= E 



(2.19) 



(2.20) 



k— — oo 



is the third Jacobi theta function ll4ll |42|] . For both distri- 
butions /i represents the main direction, while a and k are 
parameters related to the concentration [43]. 

From the previous checklist, the wrapped Gaussian satis- 
fies properties 1 to 3, while the von Mises satisfies property 4 
when the variance is replaced by its circular version (so it rep- 
resents the minimally prejudiced angle distribution, given the 
information constraints |44]). Therefore, it is tempting to side 
with the former Additionally, the Jacobi -d^ function is the 
solution of the diffusion equation on a circle with the initial 
state being a delta function, which is another way of defining 
a Gaussian wave function ifisll . However, note that if we take 
such a route, Gaussian wave functions do not lead to Gaussian 
probability distributions anymore (because the square of a 1)3 
is not a 1^3 function), a limitation that does not apply to the 
von Mises. 

In this respect, it is convenient to make a small detour into 
the question of coherent states ['46^ (recall that for the har- 
monic oscillator they are precisely Gaussian wave packets). 
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Possible definitions of coherent states for a particle on a cir- 
cle have been outlined in the literature ll47i l48ll . but they are 
of very mathematical nature. We prefer to adopt the ideas 
of Rembieliriski and coworkers [49] and construct coherent 
states \w) as eigenstates of the operator 



1/2 



E, 



so that 



W^jw) = w\w) , 



(2.21) 



(2.22) 



where the complex number w = e parametrizes the unit 
cylinder. Note in passing that 



W\m) = e™-i/2|m-l), 



(2.23) 



packets for the variables (12.261 1. In fact, they are also mini- 
mum for more intricate uncertainties [23]. However, the prob- 
lem is that, at difference of the harmonic oscillator, we do not 
have any clear operational prescription of how to measure the 
quadratures (12.26b . so they give no real physical insight into 
the statistical description of angle. 

Let us then turn back to the general commutation relation 
(12.4b . First, we observe that dealing with angle mean and vari- 
ance in the ordinary way has drawbacks. Consider, for exam- 
ple, a sharp angle distribution localized at the origin and the 
same one shifted by tt. Despite the fact that the physical infor- 
mation they convey is the same, in the later case the variance is 
much bigger Since angle is periodic but variance is not, it has 
little meaning to consider the angle measurement itself [50]. 

In circular statistics one usually calculates the moments of 
the exponential of the angle |51, ^2,^3j^], that are referred 
to as circular moments and give rise, e. g., to a circular vari- 
ance 



with [W, W''] = sinh(l)e^^. The projection of the vector \w) 
onto the basis |to) gives then 



W 



(2.24) 



while in the angular basis the corresponding expression is 



1 



(2.25) 



where ii = 6 + i£. 

In consequence, we have found three families of states with 
interesting properties: (i) states with von Mises probability 
density, Eq. ( 12.17b ; (ii) states with wrapped Gaussian prob- 
ability distribution, Eq. ( 12.18b : and (iii) coherent states with 
wrapped Gaussian amplitude density, Eq. (12.25b . Neverthe- 
less, leaving aside fundamental reasons, for computational 
purposes these three families have very similar angular shapes 
and give almost indistinguishable numerical results. There- 
fore, sometimes we will use von Mises states because of their 
simphcity and the possibility of obtaining analytical results. 



C. Constrained intelligent states 

Coherent states for the harmonic oscillator are also min- 
imum uncertainty wave packets. Given the analogy of W 
in Eq. ( 12.231 ) with the standard annihilation operator, one is 
tempted to introduce quadrature-like combinations 



P=-^{W-W'^), (2.26) 
V 2i 



that satisfy the uncertainty principle 



1, 



(AQ)2(AP)2>-|([Q,P])|^ 



(2.27) 



where (AA)^ — {A^) — {A^)"^ is the standard variance. A 
lengthy calculation |3] shows that the coherent states (12. 25b 
obey ( 12.27b as an equality and so they are indeed minimum 



al = l-\{e^^)\\ 



where 



i>P\ 



(2.28) 



(2.29) 



and p{(p) is the probability density. It possesses all the good 
properties expected: it is periodic, the shifted distributions 
P{(p + (j)') are characterized by the same resolution, and for 
sharp angle distributions it coincides with the standard vari- 
ance since |(e*'^)p ~ 1 + (0^). Moreover, this circular vari- 
ance coincides with 



(2.30) 



which is the natural extension of variance for unitary opera- 
tors Ii34.] . 

If we use ( 12.30b . the uncertainty relation associated with 
(|24] | reads 

(A^)2(AL)2 > i[l-(A^)2]. (2.31) 

Sometimes it probes convenient to express this in terms of the 
corresponding Hermitian components C and S. We have 



while [C,S] = 0, so that 



(A^)2(AL)2>i|(5)|^ 



(2.32) 



(2.33) 



Both inequalities depend on the choice of state used to evalu- 
ate (C) and (S). So intelligent states need to be distinguished 
from minimum uncertainty states: there are intelligent states 
for which the right-hand side of Eq. (12. 33b is not the obvious 
minimum value of 0. The condition of intelligence for, say, 
the first of (12.33b . reads as 



(2.34) 
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FIG. 1: Plot of the functions ceo(7;, q) (top) and 062(17, q) (bottom). 
On the right, we show two-dimensional sections of these functions 
for the values q = 0, 1, 2, 3, 4, and 100. 



that once expressed in the angle representation can be imme- 
diately solved to give 



1 



exp(i/i(/) + Ksinc 



so that the associated probability is the von Mises distribution. 
The intelligent states for the the second equation in ( 12.33b can 
be worked out in the same way. However, it is not difficult 
to prove that both inequalities cannot be saturated simultane- 
ously 1 55]. In other words, the fundamental relation (12.31b 
can never hold as an equality: it is exact, but is too weak. 

To get an attainable bound we look instead at normalized 
states that minimize the uncertainty product (AE)'^ (^^)^ ci- 
ther for a given (A_E)^ or for a given (AL)^, which we call 
constrained intelligent states. We use the method of unde- 
termined multipliers, so the linear combination of variations 
leads to 



[L^+pL+{q*E + qE^)/2]\^)^a\'f), 



(2.36) 



where p, q, and a are Lagrange multipliers. Working in the 
angle representation, the change of variables exp{ip(f))'^ (cf)) 
eliminates the linear term from (12.36b . In addition, we can 
take q to be a real number, since this merely introduces a 
global phase shift. We finally get 



- 2(jcos(2?7)] *(77) = 0, 



where we have introduced the rescaled variable 77 = 0/2, 
which has a domain < 77 < 27r and plays the role of po- 
lar angle in elliptic coordinates. Equation ( 12.37b is precisely 
the standard form of the Mathieu equation, which has many 
applications not only in optics, but also in other branches of 
modern physics [56]. An uncertainty relation of this type has 
been already investigated by Opatrny |57]. In our case, the 



only acceptable Mathieu functions are those periodic with pe- 
riod of TT or 2ti. The values of a in Eq. ( 12. 37b that satisfy this 
condition are the eigenvalues. We have then two families of 
independent solutions, namely the angular Mathieu functions 
ce„(77, q) and se„(77, q) with n — 0, 1, 2, . . ., which are usu- 
ally known as the elliptic cosine and sine, respectively. The 
parity of these functions is exactly the same as their trigono- 
metric counterparts; that is, ce„(77, q) is even and se„(77, q) is 
odd in 77, while they have period tt when n is even or period 
27r when n is odd. To illustrate these behaviors, in Fig. [1] we 
have plotted wave functions cen{rj,q) of orders n = Q and 

72 = 2. 

Since the 27r periodicity in (p requires tt periodicity in 77, 
the acceptable solutions for our eigenvalue problem are the in- 
dependent Mathieu functions ce2„(77, q) and 862,1(77, q), with 
n = 0, 1, . . .. In what follows, we consider only even so- 
lutions 062,1(77,(7), although the treatment can be obviously 
extended to the odd ones with analogous results. We take then 



*2„(?7,g) 



ce2„(»7,g) , 



(2.38) 



where we have made use of the property 

C6m 

(77, q) 66,1(77, q)dr] = ttS 



(2.39) 



(2.35) normaUze the wave function. Using ( 12.381 ) we have 



{AL)l 



1 f 

2^ Jo '^^ "^'^2n('?,9) 

][Ag-\q)-2qe,,M] 



(2.40) 



(A^)L = 1- 



- / dri cclJr],q) 003(277) 



= i-|e2„(g)r 

To obtain these analytical expressions we have expanded 
ce2n('7, q) in Fourier series 

00 

ce2„(r7, q) = Yl 4f '(9) 003(2^77) , (2.41) 

fc=0 

and integrate term by term, in such a way that 

00 

k=0 



(2.42) 

The coefficients ^2'^"' determine the Fourier spectrum and 
(2.37) satisfy recuiTence relations that can be efficiently computed by 



a variety of methods [58]. In Fig.|2]we have plotted (AL)2 
and {AE)2n as functions of the Lagrange multiplier q. 



D. Intelligent states: asymptotic limits 

To better understand Fig.|2]we first concentrate on the limit 
of small q [i. e., large (Ai?)2„]. We expand 062„(77, g) in 
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ce,(Tl,q) 
ce (Tl,q) 



4000 



8000 



W 05 



cej(tl,q) 



FIG. 2: Plot of ( AL)i„ and ( AS)i„ for the three lowest order wave- 
functions ^'2n('7, q) in terms of the parameter q. 



powers of q and retain only linear terms IS 



ceo(?7,g) = — = 



1 

3 



1 - I cos(277) 



ce2{r],q) = cos(2r7) - q 



cos(47]) 1 
4 



12 



cos(2,)-|^* + 2),] 



cos[(2n - 2)r]] 
2n- 1 



2n+ 1 

n > 2. (2.43) 



This leads to 



5q^ 



25g^ 
144 



0{q'), 



(2.44) 
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(Ai?)^„ = 1- 



4(4n2 _ 1)2 



+ 0(9"), n^l 



showing a quadratic behavior which can be appreciated in 

Fig.m 

In the opposite limit of large q [small (Ai?)|„lwe take the 
approximation in terms of Hermite polynomials 115811 

ce2„(77, q) (X e~^"'^H2n + 0{q-^'^) , (2.45) 



where u = Iq^l^ cos r\. Apart from constant factors, the states 
( 12.451 1 look like the harmonic oscillator wave functions and we 
can use them to evaluate analytically the variances. The final 
expressions involve the modified Bessel functions /fc ( ^Jq), but 
the crucial fact is that the following simple asymptotic expres- 
sions hold: 



(2.46) 



+ Oiq-') 



showing a square-root behavior that is also apparent from 
Fig. |2] The range of moderate values of q, where the tran- 
sition between the quadratic (small q) and the square-root 
(large q) regions happens, is magnified in the inset. Accord- 
ing to Eq. (12.461 1. with increasing q the uncertainty product 
(AL)2n(A£')2ri approaches a constant value depending ex- 
clusively on the mode index n; lim^_^oo(Ai)2n(A_E)2ri = 
(4n + l)/2. These asymptotic limits, confirmed in Fig. |2] 
identify the fundamental mode n = as the minimum uncer- 
tainty state for all the values of the parameter q. Henceforth, 
we always refer to the fundamental Mathieu mode, unless the 
mode index is explicitly given. 

Finally, let us note that, from Eq. (12.43b . it follows immedi- 
ately that when q ^ the probability distribution for this fun- 
damental mode is oc [l — qcos{(j))/2\'^ ~ exp(— (7C0S 0); 
while when q ^ 00, according to Eq. (12.451 1. we have 
Hoiu/V2) = 1 and cx exp[-2^cos2(0/2)] oc 

exp(— ^cos (j>). We therefore get the interesting result that 



p((/>) cx |ceo(?7,g)|^ 



0, 



p-\/9cos(/) 



(2.47) 



and hence the optimal states with very sharp and nearly flat 
angular profiles attain the von Mises shape. 



E. Suboptimal states 

Up to now we have investigated extremal states that will be 
used in the experiments as an ultimate calibration to asses the 
performance of our setup. Here, we compare these extremal 
states with suboptimal ones. There are a plenty of possible 
candidates for that: we will select a few examples that can 
be easily prepared and intuitively grasp various features of "a 
well localized angle". 

The wedge structure is our first representative. The aper- 
ture function possesses sharp edges and may be defined in the 
angle representation as 




101 < a/2. 
|(/)| >a/2^ 



(2.48) 



a being the opening angle of the wedge. The probability dis- 
tribution of angular momentum pm = l^'mp can be calcu- 
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lated using Eq. ( 12.71 ) and one finds 

(AL)2^oo, 

a sin^(mQ/2) 
^"=2^ W2)2 ' [^^'^g^] (2.49) 

(A£;)2 = l-4sin2(a/2). 

Similarly to the Fraunhofer diffraction pattern observed be- 
hind a rectangular slit, the variance ( AL)-^ is infinite due to the 
heavy tails of the sine distribution and thus the angle-angular 
momentum uncertainty relation is trivially satisfied for ( 12.481 ). 
In spite of this divergence, the experimenter cannot establish 
this simple fact from a finite data set: in general, the sampled 
angular-momentum uncertainty will grow with the size of data 
acquired. Most of the detected m fall within the central peak 
|m| < 27r/a of the distribution, which tends to regularize the 
unbounded uncertainty product. 

As another candidate for a simple single-peaked angular 
distribution we choose the state 




0.4 0.6 0.8 0.8 0.9 

AE AE 



FIG. 3: Theoretical uncertainty products for angle and angular- 
momentum variables calculated for (a) wedge angle distribution, (b) 
cosine distribution; (c) von Mises distribution, and (d) truncated nor- 
mal distribution. In all panels the solid line denotes the optimal un- 
certainty product generated by the intelligent Mathieu wave function. 
In panel (a) the broken and dotted lines correspond to pm truncated 
at the first and second minimum, respectively. 




— cos(0/a), 
7ra 



|(/.| <7ra/2, 
101 > 7rQ;/2, 



(2.50) 



that is, the positive cosine half-wave stretched to the inter- 
val — 7ra/2 < < 7ra/2, so that a < 2. The following 
discussion is valid even for a > 2 provided the cosine half- 
wave, which now spans an interval of width larger than 27r, 
is wrapped onto the unit circle. Since the delimiting aperture 
has no sharp edges, one can expect more regular results. A 
straightforward calculation yields 



(AL)2 = 1/c 



(Ai?)2 



4a cos'^(7rmQ;/2) 

^ 64sin2(7rQ!/2) 
7r2Q;2(Q;2 - 4)2 



[cosine] 



(2.51) 



At first sight, the angular- momentum distribution in Eq. (12.51b 
has a sinc-like shape with infinitely many side lobes, which 
strongly resembles the one in Eq. ( |2.49t . However, here 
the higher-order contributions to the angular-momentum vari- 
ance are negligible and both uncertainties appear to be fi- 
nite. Notice that the parameter a has now a very simple 
physical meaning: it is inversely proportional to the angular- 
momentum uncertainty. 

Next, we consider the von Mises wave function 



1 



„cos0/(2q) 



(2.52) 



v/2^/o(l/a) 

a being a monotonic function of the angular width. According 



to Eq. ( 12.471 ). the relevant uncertainties, 

h{\/a) 



Pm 



4a/o(l/a) ' 
/o(l/«) ' 



[von Mises] 



(2.53) 



1 



can be also taken (with an appropriate fitting of the parame- 
ter a) as excellent approximations to the uncertainties of the 
optimal Mathieu states, especially in the regions of small and 
large variances. 

Finally, we evaluate the uncertainties of the truncated Gaus- 
sian 



(2.54) 



which minimizes the uncertainty product for angular mo- 
mentum and angle, the latter in the sense of ordinary vari- 
ance 1I25I I26I1 . Such states become suboptimal when the vari- 
ance is replaced by a periodic measure, such as the circular 
variance advocated in this paper: 



(AL)^ 



a 
Y 



1 

2 / 

g-m /o 



erf (vra) 



"^{Rc[erf(^)]} 



(AE)^ = 1 



[truncated] 



erf (ira) 



(2.55) 



where erf (z) is the error function. 
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Figure [3] shows the comparison of these four states with 
the optimal Mathieu one. As it has been already mentioned, 
the measured uncertainty product for the wedge distribution 
grows with the size of the data acquired, since more and more 
side maxima are sampled [see panel (a)]. On the other hand, 
the uncertainty product of the cosine distribution in panel (b) 
is well defined and lies well above the quantum limit given 
by the Mathieu profile. In agreement with the asymptotic 
analysis of the previous section, the uncertainty product for 
the von Mises angular distribution falls very close to the op- 
timal curve: only at intermediate angular spreads AE we see 
a significant deviation from the standard quantum limit, while 
for the truncated Gaussian states the deviation is larger and 
shifts toward higher values of AE. The point is whether the 
currently available measurements have sufficient resolution to 
discriminate between the optimal and suboptimal states men- 
tioned above. This question is addressed in the next two sec- 
tions. 



III. MEASUREMENT OF THE ORBITAL ANGULAR 
MOMENTUM 

A. Single vortex beam 

The theory presented thus far can be applied to a variety of 
physical systems. Here, we consider a particularly appealing 
realization of a planar rotator in terms of optical beams. 

Light beams can carry angular momentum, which com- 
prises spin and orbital components that are associated with 
polarization and helical-phase fronts, respectively. In gen- 
eral, the spin and orbital contributions cannot be considered 
separately, but in the paraxial approximation both contribu- 
tions can be measured and manipulated independently. We 
emphasize that this OAM manifests at the macroscopic and 
single-photon levels and therefore paraxial quantum optics is 
the most convenient context in which to treat the OAM of light 
as a quantum resource. 

In consequence, we can leave aside the spin part and con- 
sider the simplest scalar monochromatic beam carrying OAM: 
this is precisely a vortex beam; i.e., a beam whose phase varies 
in a corkscrew-like manner along the direction of propagation. 
The corresponding spatial amplitude can be written as 



U{r) = u(r) exp{im(f>) 



(3.1) 



where we have assumed that the beam propagates dominantly 
along the z axis, so we have cylindrical symmetry. According 
to the representation in Eq. ( 12.31 ), ( 13.11 ) is an eigenstate of L 
with eigenvalue m, which is also known as the topological 
charge (or helicity) of the vortex. To check this interpretation, 
note that, the OAM density is also dominantly along the z 
axis and is given by / = rS^ jc? , where is the azimuthal 
component of the Poynting vector In a scalar theory, the time- 
averaged Poynting vector can be calculated computed as S = 
isoUj{U*WU — UWU*), where sq is a constant (with units of 
m s). The density of the OAM of the vortex beam ( 13. It then 
depends on its intensity / = |C/p and wave front helicity and 



can be expressed in terms of its power P as 

2LusomP 
1^ = — . 



(3.2) 



If we divide now by the total energy density of the field we 
finally get that the OAM per photon can be interpreted pre- 
cisely as the topological charge m. In this way, light beams 
prepared in OAM eigenstates can be used in quantum optics 
experiments in the same way as qudits. 



B. Principle of the measuring method 

A general scheme of our experimental method is sketched 
in Fig. m A collimated Gaussian beam with complex am- 
plitude Ug illuminates an amplitude mask (with transmis- 
sion coefficient Ia) performing an angular limitation of the 
beam. Immediately behind the mask, the beam transverse 
profile has a cake-slice shape given by Ua = ^aUg- Ac- 
cording to Eq. (I2.7l i. the field azimuthal amplitude distribution 
results in a spread of the spectrum composed of vortex compo- 
nents with different topological charges and amplitudes. The 
beam propagates toward a spiral phase mask [with transmis- 
sion tp = cxp{—iN(j>)] introducing a helicity —N. Behind 
the phase mask, the transmitted field is Fourier analyzed, so 
its spatial spectrum Up can be obtained as 



Up = TiitAUG * h)tp] , 



(3.3) 



where JF denotes the Fourier transform, h is the impulse re- 
sponse function of free-space propagation between the ampli- 
tude and phase masks and * is the convolution product. The 
detected power can be thus interpreted as the spectral intensity 
collected at the power meter placed at the back focal plane of 
the lens used for the optical implementation of the Fourier 
transform. 

As we show below, provided the aperture radius of the 
power meter is suitably chosen, the measured power can be 
used for estimating the OAM of the vortex mode with topo- 
logical charge N in the field behind the amplitude mask. 
This procedure can be repeated with spiral phase masks of 
different topological charges yielding the vortex distribution 
(i.e., angular-momentum spectrum) of the prepared angular- 
restricted beam. Several experimental realizations of this idea 
have been proposed and realized ll26l l27i l28ll . differing in 
technical details and data analysis [o9l l6C)ll . 

To put this in quantitative terms, let us first introduce a 
mode decomposition of the amplitude mask 



tA = 



m— — oo 



, exp(zm(/ 



(3.4) 



where a™ are Fourier coefficients. We assume that the waist 
of the Gaussian beam (of width wq) is placed exactly at the 
mask plane. The transmitted field propagates through free 
space, so the complex amplitude Uj of the field impinging 
on the spiral phase mask can be written in the form 



Uj{r) = 



amUjnir, z) e-Kp{imq 



(3.5) 
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FIG. 4: Principle of tiie measurement of the angular momentum spectrum of the angular restricted field. 



where we have used cylindrical coordinates (r, (/), z) and 

Um{r, z) = 27ruo«™^o exp ^^(r, z) , 

A,r,{r,z)= exp{-ar'^)JUPrr'ydr' , (3.6) 
Jo 

i 1 ik 3 ^ 

ho = — , = — + TT ^ P = - ■ 

AZ Wq ZZ Z 

Here Jm and uq denote the m-th order Bessel function of the 
first kind and a constant amplitude of the Gaussian beam, re- 
spectively, and z is the distance between amplitude and phase 
masks. The integration in (13.61 ) can be carried out and results 
in 

A,n{r,z) = r^,/-exp(-(5r2) 
p \ a 

X [/i(,„_i)(gr2)-/i(™H.i)(Qr2)], (3.7) 

where Q = (3^ /{8a). 

After transmission through the spiral phase mask, the 
Fourier transform of the field is performed optically and the 
spatial distribution at the back focal plane of the Fourier lens 
can be represented by the complex amplitude U f given by 

OO 

+ X! °'rnVm{v) exp[i(m - N)'4>\ , (3.8) 



m— — OO 



where 



2tt / Umir, z)Jo{2T:i'r)rdr , 



(3.9) 



Vm{i^) = 27ri'^'" / Um{r,z)Jm-N{2-Kvr)rdr . 



Here {v, tp) are polar coordinates in the transverse Fourier 
plane {xftUf), defined as = \/ x\ + yp/{Xf') and ip ~ 
a.rcta.n{yF /xp), /' being the lens focal length. The power 



captured by the circular aperture of the meter placed at the 
focal plane of the lens is given by 



P 



Jo 



\UF{v,i^)\^vdiyd^ , 



(3.10) 



where i^o = R/ (A/') and R denotes the aperture radius. Sub- 
stituting ( 13.8b into (13.10b . the detected power can be expressed 
as 



where 



P^Pn+Pc, 



Pn = 27r|aArp / \uN\'^i'dv, 



(3.11) 



Pc^ 271 ^|a,„|M 

tin rv^ 



\v.m\ Vdv . 



(3.12) 



va—~oo 



This power appears as composed of two terms: the first term, 
P/v, represents the power carried by the vortex mode of topo- 
logical charge N . The second term, Pc, is a crosstalk that rep- 
resents disturbing contributions of the remaining vortex com- 
ponents and it can be reduced by a convenient choice of the 
aperture radius R of the power meter. This possibility fol- 
lows directly from Eq. ( 13.8b . since |u7v(0)p ^ and has 
its maximum in the middle of the receiving aperture, while 
|i^m(0)P — for m ^ N . The spectral intensity of the vortex 
whose power is measured creates a sharply peaked bright spot 
at the Fourier plane, meanwhile, the spectral intensities con- 
tributing to the crosstalk power have an annular form. These 
different spatial shapes enable an optimal choice of the radius 
of the receiving aperture. In this case, both the power lost 
of the measured vortex and the influence of the crosstalk are 
minimized. 



IV. EXPERIMENTAL RESULTS 

To verify our theory, the angle-angular momentum uncer- 
tainty products were experimentally measured on various light 
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FIG. 5: Experimental setup for the generation of beams with an ar- 
bitrary transverse profile and subsequent detection of their angular 
momentum spectrum. 



beams. Given the small difference between the optimal Math- 
ieu beams and other suboptimal single-peaked angular distri- 
butions, such a measurement is also an indicator of the resolu- 
tion attainable with the present commercially available tech- 
nology. 

Figure|5]shows our setup. The beam generated by the laser 
(Verdi V2: 532 nm, 20 mW) is spatially filtered, expanded 
and collimated by the lens and impinges on the hologram gen- 
erated by the amplitude spatial light modulator (SLM) (CRL 
Opto, 1024 X 768 pixels). The bitmap of the hologram is com- 
puted as an interference pattern of the tested state (with the de- 
sired angular amplitude distribution) and an inclined reference 
plane wave. After illuminating the hologram with the colli- 
mated beam, the Fourier spectrum of the transmitted beam is 
localized at the back focal plane of the first Fourier lens FLi. 
It consists of three diffraction orders (—1, 0, +1). The unde- 
sired and — 1 orders are removed by a spatial filter. After in- 
verse Fourier transformation, performed by the second Fourier 
lens FL2, a collimated beam with the required complex ampU- 
tude profile Ua = tAUc is obtained. This completes the state 
preparation. 

The analysis begins by reflecting the prepared field Ua at a 
phase SLM (Boulder, 512 x 512 pixels), whose reflectivity is 
proportional to tp (x e~'^^'^. As it has been discussed in the 
previous section, after the Fourier transformation of the re- 
flected field, the spectral component whose helicity was elim- 
inated by the phase SLM gives rise to a bright spot (Fig.|5t), 
while the other components have an annular intensity distri- 
bution (Fig. The vortex components of the spiral spec- 
trum can be subsequently selected by the phase SLM and 
their 0AM determined by a power measurement performed 
with an optimal aperture size of the power meter To suppress 
crosstalks, the calibrating response functions were acquired 
for each phase mask. 

After the setup was carefully aligned using Laguerre-Gauss 
beams, transverse amplitude distributions of different shapes 
and angular variances were generated. Each beam was then 
scanned for values of helicities in the range of m e [—15, 15]. 
A typical transverse intensity profile and the corresponding 
measured raw data are shown in Fig.|6] 

In addition, the response functions were measured for pure 
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FIG. 6: Preparation and measurement of Mathieu states having cir- 
cular variances (from top to bottom) AE = 0.31,0.54,0.79, and 
0.91. Left: computed intensity distribution; right: measured angular- 
momentum spectrum (in arbitrary units). 
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FIG. 7: Measured response functions of our detection scheme for 
pure vortex modes of helicities m = 0, 1, 2, and 5. 



vortex modes (see Fig. |2ll. For an ideal detection, the pure 
vortex mode with the topological charge should have a S- 
like angular momentum spectrum with a sharp peak at to = 
A^. Any real detection scheme suffers from crosstalks be- 
tween modes, which tends to broaden the measured spectra 
in Fig.|6j this effect becomes more pronounced for larger he- 
licities (m > 4). It can be seen, from comparing Figs. |6] and 
|7] that the reliability of the measured spectrum decreases from 
the center (to — 0) to the edges and that beams of smaller vari- 
ances have broader angular spectra and vice versa. Hence the 
reliability of experimentally determined uncertainty products 
is expected to increase with variance. 
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FIG. 8: Analysis of measured angular momentum spectrum of a 
Mathieu beam of AE — 0.79. Raw data (small black circles), de- 
convoluted data (large gray circles) and best fit with a theoretical 
distribution (+ symbols) are shown. 



In the next step, the acquired response functions were used 
to increase the resolution of our detection scheme. Since 
the measured spectra could be considered to be convolutions 
of true angular spectra and known response functions, we 
could apply an inverse transformation to minimize the effect 
of crosstalks. Then, the angular-momentum variances were 
estimated by fitting the deconvoluted data to the theoretically- 
calculated distributions. The family of distributions used for 
fitting experimental data was parametrized by an overall nor- 
malization factor and a parameter characterizing the angular 
width of the corresponding state. For example, the fitting pro- 
cedure applied to data measured on a Mathieu beam yielded 
the value of parameter q, which was then used to determine 
the variance of angular momentum via Eq. ( I2.41l i. This stage 
of analysis is illustrated in Fig. [8] Besides getting (AL^) the 
quality of the best fit was quantified enabling to place error- 
bars on the resulting uncertainty products. 

Experimental results are summarized in Fig. |9] Given the 
resolution of the setup (indicated by error bars), the obtained 
uncertainty products fit quite well the theoretical predictions. 
As anticipated, the resolution is not uniform and gets better in 
the region of large variances. 

The inspection of the upper panels of Fig. |9] shows that the 
measurements of the optimal Mathieu and von Mises beams 
yield very similar results. This could be expected, since 
the difference between the uncertainty products of these two 
beams (see Fig. [3};) is below the resolution of the present 
setup. The cosine and wedge angle distribution in the bottom 
panels of Fig. |9]can be discriminated from the optimum more 
easily. While the suboptimality of the cosine distribution is 
confirmed only for moderate to large variances AE > 0.75, 
the wedge angular shape shows entirely different behavior: 
the uncertainty product increases with the variance. This ten- 
dency of the wedge distribution can be readily explained: as 
the variance gets larger, more and more side maxima of the 
sinc-like angular momentum spectrum fall into the detected 
window me [—15,15], yielding a corresponding grow of the 
uncertainty product. 
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FIG. 9: Experimentally determined uncertainty products for angle 
and angular momentum. The following angular distributions of beam 
amplitude in the transverse plane were measured: (a) Mathieu dis- 
tribution of Eq. ( 12.38b . (b) von Mises distribution of Eq. l l2.52t . 
(c) cosine distribution of Eq. ( I2.50t . and (d) wedge distribution of 
Eq. J2.48t . Experimentally obtained uncertainty products are de- 
noted by circles. For comparison, theoretical uncertainty products 
of the optimal Mathieu angular distribution (solid line), cosine angu- 
lar distribution (broken line), and wedge angular distributions whose 
angular momentum spectra have been truncated at the first, second, 
etc. minimum (dotted lines) are also shown. 



V. CONCLUSIONS 



In conclusion, we have formulated rigorous uncertainty re- 
lations for angle and angular momentum based on circular 
variance as a proper statistical measure of angular error. Fun- 
damental Mathieu states were identified as intelligent states 
under the constraint of given uncertainty either in angle or 
in angular momentum. In this sense, the Mathieu states pro- 
vide the optimal distribution of information between the two 
observables with possible applications in information process- 
ing. An optical test of the uncertainty relations was performed 
by using spatial light modulators both for the beam prepara- 
tion and analysis. 

Although the present experiment nicely confirmed our the- 
ory, the resolution of the present setup was not sufficient for 
observing finer details in the angular-momentum representa- 
tion of light beams. Further improvements both on the de- 
tection scheme and on hardware are highly desirable. Our 
scheme is conceptually simple but suffers from crosstalks and 
artifacts, especially at large helicities. New detection schemes 
based on direct sensing of beam wavefront could perhaps 
solve this problem. Concerning beam manipulation, spatial 
light modulators used in our experiment, though very flexible 
and easy-to-use devices, have also their drawbacks, namely 
small light efficiencies and pixellated structures. A possible 
future upgrade of the experimental setup lies in employing the 
optically-addressed SLM. 
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